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Abstract 

The Yano-Koonin-Podgoretskii (YKP) parametrisation of Hanbury 
Brown-Twiss (HBT) two-particle correlation functions opens new strategies 
for extracting the emission duration and testing the longitudinal expansion 
in heavy-ion collisions. Based on the recently derived model-independent ex- 
pressions, we present a detailed parameter study of the YKP parameters for a 
finite, hydro dynamically expanding source model of heavy-ion collisions. For 
the class of models studied here, we show that the three YKP radius parame- 
ters have an interpretation as longitudinal extension, transverse extension and 
emission duration of the source in the YKP frame. This frame is specified by 
the fourth fit parameter, the Yano-Koonin velocity which describes to a good 
approximation the velocity of the fiuid element with highest emissivity and 
allows to test for the longitudinal expansion of the source. Deviations from 
this interpretation of the YKP parameters are discussed quantitatively. 

PACS numbers: 25.75.Gz,25.75.Ld,12.38.Mh 
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I. INTRODUCTION 



The spatio-temporal extension and evolution of the interaction region in heavy-ion colli- 
sions are not directly observable. Indirect experimental access to its geometry and dynamics 
is possible through Hanbury Brown-Twiss (HBT) intensity interferometry However, 
the interpretation of the measured HBT correlations is in general model dependent, and the 
question arises to what extent their interpretational ambiguity can be reduced by a refined 
analysis of the data. 

In general, HBT radius parameters measure the Gaussian widths (second central mo- 
ments) of the source distribution in space-time [§-0] • It is the finite lifetime of the particle 
source in heavy ion collisions which complicates their interpretation. For a static boson 
emitting source, the HBT-radii have a unique interpretation in terms of geometrical source 
sizes. For dynamical sources like those created in heavy ion collisions, however, the HBT- 
radii measure certain linear combinations of the lifetime, the geometrical sizes and other 
space-time correlations Furthermore, if the source expands all HBT parameters 

become functions of the pair momentum P,^. 

For azimuthally symmetric sources, corresponding to heavy-ion collisions at zero impact 
parameter, there exist two different "complete" Gaussian parametrizations for the correlation 
function: the Cartesian parametrization with parameters Rg, Ro, Ri and Roi and 
the Yano-Koonin-Podgoretskii (YKP) parametrisation with parameters R±, R\\, Rq and 



10|-[12|,p|,|7|| . In each case, for expanding sources, these parameters are functions of the 



pair momentum. The i?-parameters have the dimension of a length while f is a velocity. 
The Cartesian parametrization has the additional difficulty that the value and the spatio- 
temporal interpretation of its parameters depend strongly on the longitudinal rest frame 
of the observer. The YKP radius parameters R±, R\\, and Rq, on the other hand, are 
independent of the longitudinal velocity of the frame in which the particle momenta are 
measured. The fourth YKP fit parameter, the Yano-Koonin velocity v, singles out a specific 
longitudinal rest frame (relative to the observer) in which the spatio-temporal interpretation 
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of the (longitudinally boost-invariant) YKP-radius parameters becomes particularly simple. 
In fact, for "transparent" sources (i.e. sources, for which particle emission occurs from the 
whole volume and is not surface dominated) without collective transverse expansion, it was 
shown in ||^ that R±, R\\, and Rq give exactly the transverse, longitudinal and temporal 
widths respectively, of the source emission function in the Yano-Koonin frame where v = 0. 
We will show here that in this case, up to small corrections from asymmetries of the source 
rapidity profile due to longitudinal expansion flow, v coincides with the longitudinal velocity 
of the fluid element around the point of highest emissivity in the source, such that R±, R\\ 
and -Ro measure (and cleanly separate) the transverse, longitudinal and temporal lengths of 
homogeneity of the source in the rest frame of the emitter. As discussed in Sec. A\ , the 



pair momentum dependence of v allows to measure the longitudinal expansion of the source 
in a very direct way. 

For the class of models studied here, we will also show (see Sec. |IV D| ) that in the absence 
of transverse collective expansion the YKP parameters show perfect M^-scaling, i.e., they 



are, for a given source, universal functions of the transverse mass M_|_ = ym^ + Kj^ of 
the particle pair, independent of the particle rest mass. This can be tested by comparing 
TTTT, KK, and pp correlations. The same is not true for the Cartesian parameters which 
contain additional kinematic and frame dependent factors which distinguish between pairs of 
particles with different mass. This type of Mj_ scaling of the YKP parameters is violated for 
sources with transverse collective expansion; also other types of transverse x-p-correlations, 



like e.g. those occurring in opaque sources [|T^, and final state effects like resonance decays 



after freeze-out [Q, can break this scaling. For a detailed discussion of these effects see 



Refs. (opaque sources) and [p!6| , p!7| (resonance decays). 

One of the purposes of this paper is to investigate the specific effects on the correlator 
and the related corrections to the spatio-temporal interpretation of the YKP parameters 
introduced by transverse expansion flow in the source. These will in general depend on the 
particular source model, and such an investigation thus must necessarily involve an extensive 
model study. We investigate here numerically a simple parametrization of the source which 



implements the finite longitudinal, transverse and temporal extension and the longitudinal 
expansion of realistic heavy-ion generated sources and contains the transverse expansion flow 
as a tunable parameter. We will see that strong transverse flow affects the interpretation 
of the YKP parameter Rq as the lifetime of the source (in its own rest frame) and thus 
renders the extraction of the duration of particle emission from the correlation data more 
difficult and less quantitatively reliable. On the other hand, it also spoils the M^-scaling 
of the YKP radius parameters which opens the possibility for an independent estimate of 
the transverse flow velocity from comparison between pion and kaon correlations. Clearly, 
such an estimate will remain somewhat model-dependent, but the mechanisms isolated in 
the present paper should still be very useful for qualitative consistency checks between data 
and theoretical interpretation. 

The investigation reported here has two aspects: an analytical and a numerical one. On 
the analytical level, we study the connection between the YKP parameters and the second 
space-time moments of the emission function, discussing the dependence of the latter on 
various geometric and dynamical features of particle emission. This discussion is largely 
model-independent; it is much more detailed than the short account given in ^ and should 
thus serve as a general basis of understanding which can be used to qualitatively anticipate 
the behaviour of the YKP parameters also for other source models than the one studied 
here. The numerical side of our study is, of course, model-dependent and our quantitative 
results must therefore be regarded with the necessary caution. 

Our paper is organized as follows: In Sec. |I|, we shortly compare the Cartesian and 
YKP parametrizations, thereby setting up our notation. In Sec. |nT| , we introduce a class of 
hydrodynamical models for the emission function. Sec. |^ contains a general discussion and 
a detailed numerical study of these models. We focus in particular on the effects of collective 
expansion flow on the YKP parameters. The main results are summarized in Sec. |V[ 
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II. HBT FORMALISM 



We shortly recall the basic relations between the emission function S{x, K), the measured 
two-particle correlation function C (q, K) , and the different Gaussian parametrisations of 
this correlator in terms of Cartesian or YKP radius parameters. We start from the relation 
||r8|,P]J^JT^ (here written down for bosons) 



Here, the emission function S{x,p) is the (Wigner) phase space density of the boson emitting 
sources [|T8i|9|,p!9| and denotes the probability that a boson with momentum p is emitted 



from the space time point x. It specifies the one-particle momentum spectrum -Pi(p) = 
EpdN/d^p = J d'^x S{x,p) as well as the two-particle correlation C(q, K). The r.h.s. of 
has to be evaluated a.t K = + P2) (the average momentum of the particle pair) and 
q = Pi — P2 (their corresponding relative momentum) where the Pi are on-shell. The Fourier 
transform in (|2.1|) does not have a unique inverse since the four components of the relative 
momentum q are not independent, due to the on-shell constraint 

n ^ ^ K K , , 

g° = /3-q, /3 = _^_, (2.2) 

which follows from q-K = 0. In practice the analysis of HBT correlation data must therefore 
be based on a comparison with specific models for the emission function S{x,K), with the 
aim of constraining the class of "reasonable" model sources as far as possible. An important 
tool for this procedure are the model-independent expressions for the HBT parameters 
which allow to calculate from an arbitrary emission function S the characteristic parameters 
of the two-particle correlation function C by simple quadrature. Experimentally, these 
HBT parameters are obtained via a multidimensional Gaussian fit to C (q, K) in momentum 
space. To compute these Gaussian parameters of the (momentum) correlation function C 
it is sufficient to use the Gaussian approximation of the (space-time) emission function S, 
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S{x, K) = N(K) S{x{K), K) exp 



6S{x,K), (2.3) 



neglecting SS{x,K) p. Here, the denote the space-time coordinates relative to the 
effective "source centre" a;(K) for pions with momentum K, 

x^iK) = x''- x'^(K), x'^(K) = {x"), (2.4) 

and 

(B-^)^,(K) = (x^x,) (2.5) 



is the inverse of the Gaussian curvature tensor in ( |2.3| ), adjusted such that the first term in 
( p.3| ) reproduces the rms width of the full source S{x^K). The (K-dependent) expectation 
values in these definitions are defined as space-time averages over the emission function: 

The correction term 5S contains information on the deviation of the emission function 
S'(x, K) from a Gaussian form in coordinate space, i.e. on sharp edges, wiggles, secondary 
peaks, etc. For the class of models discussed in this paper, however, the contributions from 
5S are known to have little influence on the half width of the correlation function [||, and 
can be neglected. Then, the two-particle correlation function C(q, K) can be calculated 
analytically from (|2.1| ): 

C(q, K) = 1 + exp [-g'^g'^(5^5,)(K)] . (2.7) 

It is fully determined by the K-dependent second space-time moments {B~^)^i, of the source 
(the "effective widths" (x^x,y)(K) or "lengths of homogeneity" [|,^]). 

A. Gaussian parametrisations of the correlation function 

In general, a Gaussian parametrisation of C(q, K) is specified by selecting a particular 
choice of three independent components of the relative momentum q and implementing in 
( [^.7D the on-shell constraint q-K = accordingly. This is usually done in a Cartesian 
coordinate system with z along the beam axis and K lying in the x-z-plane. One labels the 
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2;-component of a 3- vector by / (for longitudinal), the x-component by o (for outward) and 



the y-component by s (for side-ward). The mass-shell constraint (|2.2|) reads 



(2.8) 



where j3± = \'Kj_\/K^ ^ \'Kj_\/ Ek denotes (approximately) the velocity of the particle pair 
transverse to the beam direction, and jSi its longitudinal component. 

The standard Cartesian parametrisation of the correlation function is obtained 
by using ( |2.8| ) to eliminate from Eq. ( |2.7| ). This determines 6 Cartesian HBT radius 
parameters Rij in terms of the variances (x^a;j,)(K) of the emission function: 



C(q, K) = 1 + exp 



i,j=s,o,l 



Rl(^) = {{xi- (3it){xj - (3jt)) , i,j = s,o,l. 



(2.9) 



For an azimuthally symmetric collision region, C(q, K) is symmetric with respect to Qs 



-Qs i. Then Rl^ = R^^ = Q and 



C(q, K) = 1 + exp -i?^(K)g2 - Rl{K)ql - R^{K)qf - 2Rl,{K)q,qi 



(2.10) 



with 



Rt(K) 



^o(K) 



RU^) 



if), 

{{x - f3j){z - (3ii)) 



(2.11a) 
(2.11b) 
(2.11c) 
(2.11d) 



An alternative way of eliminating the redundant component of q in ( |2.7| ) leads to the Yano- 
Koonin-Podgoretskii parametrisation [pi [TO| , [ri| of C(q, K), 



C(q, K) = 1 + exp [-i?i(K) ql - R^K) {qf - (q^) - (r^K) + i?J(K)) (g ■ [/(K))^ 



(2.12) 



This is based on replacing in Eq. { \i.7\) qo and in terms of g_L = yg^ + Isy 'fi qi. Here, 
f/(K) is a (i^-dependent) 4-velocity with only a longitudinal spatial component: 



[/(K)=7(K)(l,0,0,t;(K)), with 7 = ^== . (2.13) 

V 1 — 

This parametrisation has the advantage that the YKP parameters R\{K), -Rg(K), and 
-R|(K) extracted from such a fit do not depend on the longitudinal velocity of the observer 
system in which the correlation function is measured; they are invariant under longitudinal 
boosts. The model-independent expressions for these YKP-parameters are most conveniently 
given in terms of the notational shorthands 

(2.14a) 

B=l{~^-^i\), (2.14b) 



/3. 



where ^ = x + iy and (y) = (xy) = for azimuthally symmetric sources such that (^^) = 
(x^ — y"^)- In terms of these expressions one finds^ 



V 


A + B 1 


~ 2C [ 




= B-vC, 


Rl 


= A~vC, 


Rl 


= m . 



2C 
A + B 



(2.15a) 
(2.15b) 
(2.15c) 
(2.15d) 



For non- vanishing transverse pair momentum Kj^, the Cartesian ( p.9| ) and the YKP ( p.l2| ) 



parametrisations are equivalent and it is instructive to compare them. The Cartesian pa- 
rameters can be calculated from the YKP ones via []^ 

Rl^ = Rl-Rl = (Rl + v^Rl) , (2.16a) 



^These expressions are valid as long as {A + B)'^ > 4C^, i.e. as long as expression ( |2.15a| ) for the 
velocity v is defined. The alternative forms of i?M and i?o given in Q are only valid if additionally 



A + B > 0. For a detailed discussion see [15|. 



Rf = {l- Pf) Rl + 7' (A - vf {Rl + i?;) , (2.16b) 
Rli = P± {-PiRl + 7' (A - v) {rI + R11)) , (2.16c) 
Rl = R\. (2.16d) 

Later we will see that, for the explicit source models studied in this paper, in most cases 
the Yano-Koonin velocity v is very close to the longitudinal pair velocity A. If this is true 
(and one should be careful not to use the following expressions without first checking this) 
Eqs. ( |2.16| b,c) simplify to 

Rl^R\h\ (2.17a) 
Rli^-(3^(3iRl. (2.17b) 



There is a slight subtlety for Kj^ = 0. In this limiting case, the on-shell constraint ( p.2|) 
reads = jSiqi and cannot be used to eliminate in (|2.7| ) go and in terms of q±, g° and 
gs. Hence, strictly speaking, the YKP parametrisation exists only for K± 7^ 0. In practice, 
however, this does not lead to complications since the Kj_ —>■ limit is well-defined for all 
YKP-parameters (see Sec. [IV B| ). 

B. Advantages and drawbacks of different Gaussian parametrizations 



The relations (|2.16|) provide a powerful consistency check on the experimental fitting 



procedure of the correlation radii. They show that both parametrizations contain exactly 
the same spatio-temporal information. However, certain space-time characteristics of the 
source are more directly accessible in a particular parametrization. This is especially the 
case for the duration of the particle emission process, the "lifetime" of the source. 



To see this we return to the expressions ( |2.11| ) for the Cartesian HBT radii. These 
mix spatial and temporal information on the source in a non-trivial and frame-dependent 
way. Their interpretation in various reference systems was analysed analytically for 
a large class of (azimuthally symmetric) model emission functions and compared with the 
numerically calculated correlation function ^j. For these models, the difference 
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Rl^ ^Rl-Rl = - 2/3^ + ((i^) - {y')) (2.18) 

is dominated by the first term on the r.h.s. and thus provides access to the hfetime At = 



(t^) — of the source Note, however, that the definitions ( |2.11| ) and ( |2.18| ) are not 



Lorentz invariant, and that the hfetime At extracted from Eq. (|2.18|) thus depends on the 



analysis frame. Furthermore, in practice the term f3j_{P) turned out to be much smaller than 
the terms (x^) and (y^) which are the leading contributions to Rl and R^, respectively 
As a consequence, excellent statistics of the data with very small statistical errors of Ro and 
Rs are required to extract the small contribution -R^jg. This makes the extraction of a small 
source lifetime from the standard fit difficult^. Successful attempts have been reported from 
low-energy heavy-ion collisions (using 2-proton correlations) where the measured lifetimes 
are very long: 25 ± 15 fm/c in Ar+Sc collisions at E/A = 80 MeV p3[ and 1400 ± 300 fm/c 
in Xe+Al collisions at E/A = 31 MeV (the latter is the typical evaporation time of a 
compound nucleus). Two-pion correlations at ultra-relativistic energies {E/A = 200 GeV) 
so far failed to yield positive evidence for a non- vanishing emission duration p5|j26|| , except 
for the heaviest collision system Pb+Pb |2^, but even there the effective lifetime is only a 
few fm / c. 

For the YKP parametrization the situation is different. The parameters Rq, R\\ and R± 
are invariant under longitudinal boosts and thus independent of the analysis frame. The key 
to their space-time interpretation is provided by the fourth fit parameter f (K). It specifies 
a (pair momentum dependent) longitudinal reference frame, the Yano-Koonin (YK) frame 
which defined by f = resp. C = (see ( p.l5a| )), in which the space-time variances (2.15) 
for the YKP radius parameters simplify considerably. Especially, for certain classes of source 
models including the one studied below, the terms proportional to {zx), {xi), and (x^ — y"^) 



^The situation may be better for very long-lived sources which are predicted by hydrodynamics 
if there is a phase transition to a quark-gluon plasma and the collision fireball is initiated within 



a certain range of energy densities [22|. 
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are small Neglecting these terms one obtains ^,0] in the YK frame 

Rl{K) = {f), (2.19a) 



2. 

1 1 



i?J(K) = 5 = ^5 - j ^ - ^(y^) ^ {z^) , (2.19b) 

2\ 



i?^(K) = A=(^(^t--5j j-^{f)^{P). (2.19c) 

In neighbouring frames, the Yano-Koonin velocity can be calculated with the same approx- 
imations as 

C (zi) 

V ^ ^ ~ ^ ' . (2.20) 

Note that in the YK frame the temporal structure of the source enters only in the parameter 



i?o- Its leading conribution is given by the time At(K) = y (t^) during which particles of 
momentum K are emitted in this frame. In the YKP parametrisation i?o ~ is fitted 
directly and not obtained as the difference of two large fit parameters as in the Cartesian 
fit. 

In practice, however, the extraction of Rq from the YKP fit is still not easy. From 
Eq. ( |2.12D it follows that in the YK frame i?o must be extracted from the g^-dependence of 



the correlator. Due to the mass-shell constraint ( |2.2D the interesting range of (f is limited, 
especially for low-momentum pairs, and the sensitivity of the fit function (|2.12|) to Rq is 



weaker than to the two other radius parameters. Rq values thus tend to come out with 
larger experimental error bars. 



III. A MODEL FOR A FINITE EXPANDING SOURCE 



For our numerical study we have taken the model from Ref. |^ with the emission function 



(27r)3.,/27r(A 



r 



K ■ u{x) 
T 



exp 



(r - Tof [t] - rioY 



2(Ar)2 2i?2 2{Ar]y 



■ (3.1) 



The first term specifies the shape of the freeze-out hypersurface, the second one is a Lorentz- 
covariant Boltzmann factor encoding the assumption of local thermal equilibration super- 
imposed by collective expansion, while the last one has a purely geometrical interpretation. 
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The space-time coordinates in longitudinal and temporal directions are parametrised by the 



space-time rapidity rj = | ln[(t + — z)] and the longitudinal proper time r = — z"^- 



In the transverse direction, the radius is r = yaP^+y^. Accordingly, the measure reads 
d'^x = T dr drj r dr d(f) . The time-component of the pair momentum is set to the on-shell 
value Kq = Ek = \/m? + K^. This approximation was studied in detail in Ref. |Q where 
it was shown to be acceptable. Thus the pair momentum K can be parametrised using the 



momentum rapidity Y = |ln[(l + /?;)/(! — A)] ^ind the transverse mass M_l = ym^ + K\, 

K^" = (M^ coshy, K^, 0, sinhY). (3.2) 

We implement longitudinal and azimuthally symmetric transverse expansion of the source 
by parametrising the flow velocity in the form 

u^ix) = [ cosh rii(T,ri) cosh rit(r), —smhrjJr), -smhrit{r), sinh?7/(r, r/) coshr^^fr) . (3.3) 
\ r r / 

For the longitudinal flow rapidity we take rii{T,ri) = rj independent of r, i.e. we assume a 
Bjorken scaling profile Vi = z/t in the longitudinal direction. For the transverse fiow 



rapidity we take a linear profile of strength rjf-. 

Vt{r) = Vf (^) • (3.4) 

The scalar product in the exponent of the Boltzmann factor generates the x-i^'-correlation 
in our source. It can then be written as 

X 

K ■ u{x) = Mj_cosh(?7 — Y) coshrit{r) — K ±— sinh rit{r) , (3.5) 

r 

Please note that for non-zero transverse momentum K±, a finite transverse fiow breaks the 
azimuthal symmetry of the emission function via the second term in ( |3.5| ). For rjf = 0, 
the emission function is azimuthally symmetric for all Kj_. Also, it then has no explicit 
i^'_L-dependence, and M± is the only relevant scale. As will be discussed in Sec. [IV D| this 
gives rise to perfect M_L-scaling of the YKP radius parameters in the absence of transverse 
fiow, which is again broken for non-zero transverse fiow pOj . 
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Besides 77/ , the model parameters are the freeze-out temperature T, the transverse geo- 
metric (Gaussian) radius R, the average freeze-out proper time tq as well as the mean proper 
emission duration Ar, the centre of the source rapidity distribution r^o? and the (Gaussian) 
width of the space-time rapidity profile Ar/. A rough spatial picture of the source at various 
fixed coordinate times can be gleaned from the Figs. 1 and 2 in Ref. (although their 
source has sharp edges whereas ours is smoothed by Gaussian profiles) and from Figs. 1 and 
2 in Ref. [l^. Note that our parametrization of S{x, K) does not allow for the case of opaque 
sources where the emission is surface dominated [l^Jl^. In this case, the contribution of 
{x^ — y'^) may become negative and large which might alter the argument following 
( p.l8| ). A detailed study of such opaque sources is presented in Ref. [|1^. On the basis of a 
comparison with the preliminary data of Ref. the authors of that study conclude that 
the source created in Pb+Pb collisions at the CERN SPS is rather transparent and not 
opaque; its quahtative features are well described by the model presented here . 



We did our calculations for pions (m = m^± = 139 MeV/c^) and kaons ( m = mx± = 494 
MeV/c^). Resonance decays |]I5|TB| are not discussed here but is deferred to a separate 
publication The calculations presented here are meant to illustrate general properties 
of the YKP parameters; no attempts to compare with data will be made. 



IV. LIFETIMES AND SIZES FROM THE YKP-FIT TO THE CORRELATION 

FUNCTION 

In this section we present a quantitative study of the YKP fit-parameters. Since the 
YKP-parameter is identical to the "side" radius of the Cartesian parametrisation, = 
R^ = {y"^), its interpretation is obvious and independent of the longitudinal velocity of the 
reference frame. Hence we focus subsequently on the remaining three parameters R^, 
and V. 

Unless stated otherwise, the numerical calculations below are done with the set of source 
parameters T = 140 MeV, R = 3 fm. At] = 1.2, tq = 3 fm/c, Ar = 1 fm/c. 
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A. The Yano-Koonin velocity 



According to Eq ( 2.13 ), the YKP fit parameter f is a longitudinal velocity. In this 
subsection we give a detailed discussion of the reference frame specified by v and establish 
its relation to several other commonly used reference frames. Their definitions are: 

• CMS: The centre of mass frame of the fireball, specified by rjo = 0. 

• LCMS (Longitudinally CoMoving System |2T|): A pion (kaon) pair-dependent frame, 
specified hj f3i = Y = 0. In this frame, only the transverse velocity component of the 
pion (kaon) pair is non-vanishing. 



LSPS (Longitudinal Saddle- Point System [0): The longitudinally moving rest frame 
of the point of maximal emissivity for a given pair momentum. In general, the velocity 
of this frame depends on the momentum of the emitted particle pair. For symmetric 
sources the point of maximal emissivity ("saddle point") coincides with the "source 
centre" a;(K) defined in ( |2.4| ). In this approximation, for a source like ( |3.1|) , the LSPS 



velocity is given by the longitudinal component of 'u'^(x(K)). 

• YK (Yano-Koonin frame 0]): The frame for which the YKP velocity parameter van- 
ishes, ^^(K) = 0. Again, this frame is in general pair momentum dependent. 

These four frames are quite different in nature. The velocities (or rapidities) of the 
CMS and LCMS frames can be easily determined experimentally, the first from the peak in 
the single particle rapidity distribution, the second from the longitudinal momentum of the 
measured pion pair. However, the velocity of the LSPS is in a sense more interesting: from 
its definition it is directly related to the longitudinal expansion velocity of the source. In 
fact, longitudinal expansion of the source leads to a characteristic dependence of the LSPS 
velocity ^Lgpg on the pair rapidity. This is most easily seen by considering two extreme 
fireball models: 
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(1) If the source does not expand, all its elements move with the same velocity (ra- 
pidity), namely that of the CMS. Then there is no kinematic difference between differ- 
ent parts of the fireball, and the saddle point is /^-independent and given by the peak 
of the space-time distribution of the source. Thus the rapidity of the LSPS, defined as 
^Lsps ~ '^LSPs)/(-'- ~ "^Lsps)]' independent of the pair rapidity Y and identical to 
the rapidity of the CMS. 

This behaviour of Y^^p^ arises automatically if the source has no x — correlations, i.e. if 
the emission function factorizes, S{x,K) = F{x) G{K). Then all space-time characteristics 
of the source, including the LSPS velocity, are determined by F{x) alone and do not depend 
on K. Hence Vlsps independent of Y. Factorization is, however, not necessary for this 
behaviour: Non- vanishing x — K correlations generated, e.g., by temperature gradients don't 
induce a F-dependence of Y^^^^ as long as the source does not expand longitudinally. 

(2) If we set in Eq. (|3.1| ) Ar] oo, we recover the Bjorken model p8|] for a longitudinally 



infinite source with boost-invariant longitudinal expansion. The only //-dependence then 
comes from the thermal Boltzmann factor, and the longitudinal saddle point obviously lies 
dXr] = Y . Since for the Bjorken scaling profile rj coincides with the longitudinal fluid rapidity, 
this implies that in this case the rapidity of the LSPS is identical with the pair rapidity, 
y^sps = ^ 1 the LSPS coincides with the LCMS. 

It is obvious from this discussion that knowledge of the function Y^^^^iY) would allow 
to distinguish between these two scenarios: A non-expanding source would yield I^lsps (^) ~ 
const, while a source with boost-invariant longitudinal expansion gives Y^^^^iY) = Y. Re- 
alistic models are expected to lie in between these two extremes. 

Unfortunately, the LSPS- velocity fLSPs(K) cannot be measured. This is clear from its 
definition as the longitudinal fiow velocity evaluated at the point of maximum emissivity. 
As discussed above, this point is approximately given by the source centre x(K) which itself 
is unmeasurable: it drops out |j^ from both the single and two-particle spectra which are 
invariant under a translation of the source centre (even if K-dependent!). The only velocity 
we can measure, namely from an YKP fit of the two-particle correlator, is the Yano-Koonin 
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(YK) velocity f (K) resp. the associated rapidity Y^^ = [(1 + '^)/(l ^ '^)]- It is therefore 
very gratifying to know that for sufficiently rapidly expanding systems the two velocities are 
always closely related (although in general not identical). In particular, we will show below 
that for the class of models the function Y^^^iY) shares with 5^lsps(^) feature that 
it provides a direct signature for longitudinal expansion. 

The close relationship between the (measurable) YK-velocity and the (theoretical) LSPS- 



velocity for certain source models has been known for years. In Ref. ||Tl| a "symmetric frame" 
was introduced as the reference frame, in which the production process is symmetric in the 
beam direction. In this frame, the longitudinal extension and the lifetime of the source reach 



their extremal values [11 1. For moving, but non-expanding azimuthally symmetric sources 



Podgoretskii found in this way the parametrisation ( p.l2| ), with K-independent parameters 



R±, R\\, and Rq, and with what we call the Yano-Koonin velocity v being identical to the 
velocity of his "symmetric frame" . In this case the YK-system also coincides with the rest 
frame of the source as a whole (CMS) as well as with the LSPS. 

That the coincidence between the YK and LSPS systems is more generally valid 
for sources which are symmetric around their saddle point x(K) has been observed in 
Refs. It thus holds for any emission function in the Gaussian saddle-point approxi- 



mation, due to the symmetry of the latter. As the following paragraph will show, differences 
between the YK-velocity and the velocity of the LSPS are only due to asymmetries of the 
source around its saddle point. Although such asjTumetries usually exist for collectively ex- 
panding sources with finite geometric extension, they are generally small and can be treated 
perturbatively. Therefore, the YK-frame and the LSPS-frame are usually very close to each 
other, V ~ "^lsps- From the examples above it is clear, however, that the same is not true 
for the LCMS (i.e. the longitudinal rest frame defined by the pair rapidity Y), and that 
generally v ^ w^cms- 

Let us now discuss the difference v — v^^p^ in more detail. If it is small, so is C when 
evaluated in the LSPS-frame. From Eqs. ( [^.20D and ( |2.14c| ) we see that then in the LSPS 
frame v is given by 
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where we expanded in first order of (tx), (zx), [x^ — y^). Tlie smallness of tliese terms 
was argued in and will be checked in the following subsection. The first and second 
term of ( ^4.1| ) reflect the longitudinal asymmetry of the source around the saddle point; they 
vanish for sources with longitudinal reflection symmetry z —>■ —z. Similarly, the second and 
third term vanish unless the reflection symmetry x — * —x around the saddle point in the 
"out" -direction is broken (e.g. by transverse flow). (Asymmetry in t is needed for the first 
and third terms to become non-zero.) A non-zero value of the last term, finally, indicates 
the breaking of the "out"- "side" rotation symmetry; this can again be caused by transverse 
source gradients as e.g. transverse flow (see Eq. ( p.5| )). We conclude that the difference 
V — f^gpg is entirely due to asymmetries of the source around the saddle point. Furthermore, 
we will show in Sec. |1V B| that the last three terms in ( |4.1| ) are small for small values of the 
transverse flow rapidity rjj and/or small values of K±. (They vanish for rjf = 0.) In these 
limits the difference between v and f lsps is dominated by the longitudinal source asymmetry, 
and V is very accurately given by the (leading) first term in ( |4.1| ) , see Eq. ( |2.2CI| ) . Note that in 



our model the breaking of the longitudinal reflection symmetry is due to the non-symmetric 
rapidity profile of the emission function for Y ^ tjq. 

For a quantitative discussion we plot in Fig. |l] the Yano-Koonin rapidity Y^^ as well as 
the difference Yy^, — ^lsps functions of M±_ and Y. All rapidities are given relative to 
the CMS. One sees that for large values of K± the agreement of with I^lsps almost 
perfect. Also, in this limit both rapidities approach the value of Y, i.e. the YK and LSPS 
systems coincide with the LCMS. The reason for this is that for large K± the Boltzmann 
term in the emission function (|3.1|) becomes sharply peaked around the point x where the 
fluid velocity agrees with the pair velocity; the geometric terms in the emission function are 
much smoother and can be neglected. The relevant term is the first term in Eq. ( p.5| ), and 
thus the relevant variable is the transverse mass M±. Hence this kinematic region starts for 
kaons at smaller values of K± than for pions (see Fig. |l|). 
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For small values of M^, the difference K^k ~^lsps increases, and both begin to lag behind 
the LCMS rapidity Y. Still, the YK frame is closer to the LCMS than is the LSPS. 

The fact that the YK and LSPS systems track each other so closely implies that the 
linear rise of the YK rapidity with the pair rapidity Y reflects nothing but a similar rise of 
the LSPS rapidity with Y. As argued above, the latter is a direct indication for longitudinal 
expansion of the source. However, it should be noted that this expansion need not necessarily 
be of hydrodynamic nature. The same feature would be generated by a source consisting 
of free-streaming pions and resonances which were created at an initial proper time Tform 



through a boost-invariant production mechanism ||2^, suffering no further re-scattering. It is 
easily seen that the strict correlations between coordinates and momenta in a free-streaming 
gas again lead to a linear dependence of the "source rapidity" 1^^. on the pair rapidity, with 
M^-independent unit slope. In fact, it is possible to simulate this situation with the emission 
function ( p.l| , p.3|) by setting T and f]f to zero, i.e. by eliminating the thermal smearing of the 
momenta and the transverse collective flow. (Of course, this would also result in vanishing 
YKP radius parameters, because pions of fixed rapidity Y can come from only a single point 
in the source.) 

A linear rise of Y^k with Y (with approximately unit slope) was recently observed by 
the GIBS collaboration in Dubna [jl2| by analysing pion correlations from Mg+Mg collisions 



at 4.4 A GeV/c. They interpreted their result as evidence for rapid longitudinal expansion 
of the source. The data were averaged over the transverse pair momentum K±. The data 
sample was taken with a "central" trigger, but since ^^Mg is a rather small nucleus with a 
large surface to volume ratio it is not clear what fraction of the participating nucleons were 
stopped to become part of a thermalized fireball. The strong linear increase of with Y 
could thus also reflect to some part the free-streaming expansion of the pion sources created 
in the periphery of the nuclear reaction. Preliminary results of the NA49 collaboration at 
CERN for Pb+Pb collisions at 158 A GeV/c also show a rise of Y^^ with the pair rapidity Y 
27| . Here, however, the y-dependence of Y^^^ does not appear to be quite linear, and its slope 



is less then 1. With all due caution with regard to the still preliminary nature of these results. 



this may indicate genuine hydrodynamic longitudinal expansion at a somewhat slower rate 
than resulting from our longitudinal scaling profile. NA49 have also looked separately at a 
subsample of pairs with K± > 300 MeV/c, showing that for them the YK rapidity appears 
to rise more rapidly with Y than in the i^'j,- averaged sample, in agreement with theoretical 
expectations for a thermalized source (see. Fig. |I]a). 

The M^-dependence of '^^^ ^e understood along the same lines. Fig. |l|b shows that 
Y^^{Y) approximates Y better with increasing transverse mass. In the limit M± — oo the 
bosons can be emitted only from the source element which moves exactly with the same 
rapidity, hence Y^k In the opposite limit M± —>■ the Boltzmann factor in ( p.l|) 

becomes a smooth function of x, and the emission function is dominated by the Gaussian 
geometric terms. In this limit one thus expects the YK rapidity to approach the value 
~ Vo- The numerical results of Fig. |l]b show that the value of the pion mass is already 
large enough for the Boltzmann part of the emission function to become important. As a 
result the M^-dependence of weak in the entire range which can be covered by pions, 

and even weaker for kaons. 

B. Correction terms 

In this subsection, we study quantitatively the correction terms of Eqs. (|2.19|) and ( [4.1|) 
which may compromise the approximation ( |2.20|) for v and the simple interpretation of R\\ 
and Rq as longitudinal and temporal widths of the emission function. Since the geometric 
interpretation of R\\ and Ro refers to the YK frame, our analysis will also be done in this 
frame. 

Let us first focus on the central rapidity region Y = r]Q. Then Y^k — ^lsps ~ i-^-' 
the four reference frames listed in Sec. |1V A| coincide. Furthermore, the source is symmetric 
in the longitudinal direction and thus (xz) = (zt) = 0. The only non- vanishing corrections 
thus arise from the terms (xt) and (x^ — y^). In Fig. |^ they are plotted as a function of K± 
for different values of the scaling parameter rjf for the transverse flow. Without transverse 
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flow (i.e. for r]f = 0) tlie source is azimutfially symmetric (see Eq. (|3.5| )) wliicli implies 
{x^ — y"^) = 0. Also, the source is reflection symmetric in the out-direction, (xt) = 0. For 
non-zero transverse flow the correction terms are generally non-zero, and they grow with 
increasing rif. Note that, for fixed rjf and K±, the correction terms are considerably smaller 
for kaons than for pions. This can be understood as follows: as discussed after Eq. (^.5[), 
for rjf = the source depends only on M±, and thus even for non-zero transverse flow 
everything to zeroth order still scales with M±. We will show below that for expanding 
sources the regions of homogeneity, which effectively contribute to the correlation function, 
are generically decreasing functions of M±. This is thus also true for the correction terms. 
Since at fixed the value of M± is larger for kaons than for pions, the corresponding 
correction terms are smaller in absolute terms (although not necessarily relative to the 
leading contributions) . 

The ii'^-dependence of the correction terms at non-zero transverse flow f]f ^ can also be 
easily understood. The rise of {x^ — y"^) for increasing transverse momentum (Fig. ^) is due 
to the azimuthal symmetry breaking by the second term of Eq. ( pj.5D which increases both 
with and rjf. It agrees with the findings of Ref. ^ but, as pointed out in |T^, contradicts 



the behaviour seen by Pratt in the first of Refs. for an infinitesimally thin spherically 
expanding shell where {x^ — y"^) decreases with increasing K±. A similar behaviour is seen 
pSi in hydrodynamical simulations where freeze-out occurs along an infinitesimally thin 
freeze-out hypersurface; there also (x^ — y'^) first decreases very rapidly with increasing K^^, 
then saturates and slightly increases again without, however, ever turning positive. The 
strong decrease of {x"^ — y"^) with appears to be an artefact of the idealization of an 



infinitesimally thin expanding shell; in ||T5[ it was shown to be much less visible for opaque 
sources with a finite thickness of the emitting surface layer, returning to the here observed 
rise already for a rather modest surface thickness. 

Different from {x^ — y'^), the variance (xt) reaches an extremum and then decreases 
again for very large K± (Fig. ^). This results from an interplay between the increasing 
breaking of the x —x reflection symmetry, which tends to increase the value for (xt), and 
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a decreasing homogeneity length in space-time rapidity t] which affects the t = t cosh rj part 
of this variance. 

In the mo del- independent expressions ( |2.19| ) and ( [4 .11 ) for the YKP parameters, the 
correction terms discussed above are divided by P± and respectively. From Appendix]^ 
we know that the ratios remain finite in the limit (3± —* 0. Still, the corrections to the YKP 
parameters could become sizeable, depending on how slowly the variances (xt) and (x^ —y"^) 
vanish in this limit. In Fig. we show that the correction terms actually remain small 
even after dividing them by the appropriate powers of f]±. Thus, at least at mid-rapidity. 



the leading order approximations ( p.l9[ b,c) are seen to be generally very good. The largest 
correction comes from the difference (x^ — in Eq. ( p.l9c|) . A more detailed discussion 

of its effects on Rq will follow in the next subsection. 

We now proceed to a discussion of the correction terms for Y ^ r]Q. We define Y^^ = 
y — r^o as the rapidity of the pair in the CMS. In Fig. ^ we compare the correction terms 
for Y^^ = 3 to those for Y^^ = 0. Of course, at Y^^ = 3 the YK frame no longer coincides 
with the CMS, see Sec. |1V A| . Since the transverse variances are not affected by longitudinal 
boosts nor do they depend on Y (see Appendix 0), (x^ — y"^) does not change with the pair 
rapidity. However, at a given value of K± the transverse pair velocity f3± = K^/Ek becomes 
smaller, since Ek increases with the pair rapidity Y^^. The correction term (x^ — y'^)/(3± 
thus increases with Y^^, especially at low K^. A similar effect is seen in the plots for 
{xi)/P±. While (xt) decreases with increasing Y^^, the ratio with P± actually increases 
by about a factor 2 at small K±. Since in the forward rapidity region the source becomes 
non-symmetric under reflection z —z, (xz) and (xz) / f3± are non-zero, but small. 

We conclude this subsection with a discussion of the sensitivity of the correction terms 
to the longitudinal extension of the source, which is parametrised by Arj. We again put 
^CM — 0' so that {xz) = 0. From the explicit expressions given in Appendix ^ it is clear 
that the transverse variances (x^ — y^) are independent of Arj. Therefore only the size of 
(xt) changes. In Fig. ^ we show its if^-dependence for different values of Arj. Due to 
our assumption of freeze-out along a hyperbola of constant proper time tq (smeared by an 
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amount At), an increase of At] causes a larger effective extension of the source both in 
the z and in the time direction. This becomes especially important at small where the 
geometrical factors in the emission function ( |3.1| ) play an important role. This explains the 
relatively large A?7-dependence of (xt) at small K± and the weaker dependence at large K± 
where the source distribution is dominated by the Boltzmann term which does not depend 
on At]. 



C. Quantitative numerical study of YKP radius parameters 

In this subsection we combine the results from the previous subsection for the correction 
terms with the leading contributions in order to arrive at a quantitative understanding of 
the longitudinal and temporal YKP radius parameters R\\ and Rq, and in particular of their 
dependence on the pair momentum K. 



In Fig. 1^ we show Rq and R\\ together with their approximations y (t^), y (5^), for pion 
pairs with rapidity Y^^ = and Y^^ = 3, respectively, as functions of Kj_. For vanishing 
transverse flow both approximations are seen to be exact, in agreement with the discussion 



from the previous subsection. For non-zero transverse flow the approximation ^ y (P) 
remains exact for pairs with Y^-^ = 0. The reason is that for such pairs the YK rapidity 
relative to the CMS is zero, and thus the longitudinal velocity jSi (which multiplies the 
correction terms in ( |2.19b| )) of the pair in the YK frame vanishes. From Fig. ^ one sees. 



however, that also for forward rapidity pairs at Y^^ = 3 the correction terms stay below 
10% at all values of K±. 

The situation is not quite as good for Rq. Here one sees apparently strong differences 



between Rq and y (t^) as soon as the transverse flow is switched on. From Fig. ^ it is clear 
that the (in our case positive) correction term (5^ — y"^) is the culprit and dominates the 
difference. For a transverse flow of rjf = 0.6 as shown in Fig. ^ this term becomes (in 
the experimentally accessible range) larger than 1 (fm/c)^ and thus comparable to the 



leading term J (t^). However, the numerical results shown in this Figure actually correspond 
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to a rather extreme situation. First, the assumed transverse flow rapidity rjf = 0.6 is large; 
the heavy ion data at AGS and CERN energies for Si- and S-induced reactions seem to 



require smaller values p4|-p^. Second, the difference between Rq and J (t^) is small at low 



transverse momenta and becomes large only at large K±] in that range the leading term 



(t^) is essentially given by the source parameter At (see discussion below) which in Pb+Pb 
and Au+Au collisions ||27|j30(| (where rif may be larger than for the smaller systems analyzed 
so far) seems to be bigger than the 1 fm/c assumed here^. 



Fig. ^a shows that the effective source lifetime At = y (t^) is a strong function of the 
pair momentum K: it is largest at small rapidity Y^^ and transverse momentum K± and 
decreases with increasing Y^^ and/or K^. Its asymptotic value for large K in the CMS is, 
not unexpectedly, given by the variance (|Allj ) of the proper time distribution of our source 
( p.l|) . But why is it larger for pairs with smaller momenta K? 

From Fig. ^ it is clear that the longitudinal region of homogeneity R\\ is a decreasing 
function of the pair momentum K. The reason for this is the same as the similar decrease of 
Ri in the Cartesian fit and well understood as a consequence of the strong longitudinal 
expansion of the source. This expansion introduces a longitudinal velocity gradient, and 
the longitudinal length of homogeneity is given by the inverse of this gradient multiplied 
by a "thermal smearing factor" [Q. The latter reflects the statistical distribution of the 
particle momenta around the local source fluid velocity, and for a thermal distribution the 
spatial region over which this thermal smearing is effective decreases with increasing pair 
momentum. This causes the shrinking of the longitudinal homogeneity length with K. 

Since for different pair momenta Ro measures the source lifetime in different YK reference 



^If the source, unlike ours, is opaque, i.e. if the particle emission is strongly surface dominated, 
(x^ — y^) tends to be negative [13|,l^j3^. In this case, the deviation of Rq from J (t^) will have the 



opposite sign |1^,33|, and Rq usually even turns negative for small values of K±. For a detailed 
study we refer to [ 15 1 . 
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frames, the freeze-out "hypersurface" will in general appear to have different shapes for pairs 
with different momenta. Only in our model, where freeze-out occurs at fixed proper time tq 
(up to a Gaussian smearing with width At), is it frame- independent. It is thus generally 
unavoidable (and here, of course, true in any frame) that freeze-out at different points z 
in the source will occur at different times t in the YK frame. Since a 2;-region of size R\\ 
contributes to the correlation function, determines how large a domain of this freeze-out 
surface (and thus how large an interval of freeze-out times in the YK frame) is sampled by 
the correlator. This interval of freeze-out times combines with the intrinsic Gaussian width 
At to yield the total effective duration of particle emission. It will be largest at small pair 
momenta where the homogeneity region R\\ is biggest, and will reduce to just the variance 
of the Gaussian proper time distribution at large pair momenta where the longitudinal (and 
transverse) homogeneity regions shrink to zero. 

Another interesting feature of Fig. ^ is that at large K± both R\\ and Rq are independent 
of the pair rapidity Y. This is a consequence of our boost-invariant longitudinal velocity 
profile and need not remain true for systems with different longitudinal expansion. As argued 
before, at large M±_ the space-time shape of the source is dominated by the Boltzmann term 
and becomes insensitive to the Gaussian geometric factors. The HBT radii thus only see 
the local velocity gradients which in our case are invariant under longitudinal boosts. At 
large M±_ pion pairs with different rapidities Y thus all see the same local source structure, 
and the YKP radii become F-independent. 

We close this subsection with a discussion of the dependence of R\\ and Rq on the other 
source parameters. Since the rapidity dependence does not change qualitatively from what 
has already been discussed, we concentrate on zero rapidity pion pairs, Y^^ = 0. For the 
transverse flow we choose a non-zero, but moderate value of rjf = 0.3. 

In Fig. ^ we show the dependence on the longitudinal size of the source which is 
parametrised by Arj. One sees that at large transverse momenta neither R\\ nor Rq and 



Y (f^) are affected by the width Ar] of the Gaussian geometric factor, in line with the ar- 

guments above. At small transverse momenta, both radii increase monotonically with Arj. 
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This means that, at low K^^ R\\ becomes sensitive to the global longitudinal geometry of 
the source and no longer only reflects the local longitudinal velocity gradients. Fig. ^ is 
very interesting: for small Ar/ the longitudinal length of homogeneity is limited by the lon- 
gitudinal geometry, and Rq never has a chance to probe a large region of the proper time 
freeze-out surface. Hence At is limited to the variance of the proper time distribution T(r) 
in the source (see ( |A2|,|A11| )), independent of K±. As A?7 increases, at small K± the longitu- 



dinal length of homogeneity increases too, and Rq receives an additional contribution from 
the time variation (in the YK rest frame) along the freeze-out surface inside a longitudinal 
region of size R\\. Thus the rise of the effective source lifetime At at small is an indi- 
rect measure for the longitudinal geometric size of the source. Unfortunately, the detailed 
quantitative dependence is model-dependent. 

Fig. 1^ shows what happens when the width of the proper time distribution T{t) 
(Eq. ([A2|) ) is changed. Increasing Ar by 1 fm/c, Rq also increases by about 1 fm/c (slightly 
less at large -R'.lQ), while R\\ increases more at small K_[^ and less at large K^. The increase 
of -R|| is due to the decrease of the longitudinal velocity gradient (which for a boost-invariant 
profile is 1/r) with r. As the time distribution T(r) becomes wider, larger proper times are 
probed by the emitted pions resulting in larger longitudinal homogeneity regions. 

Changing the average freeze-out time tq rather than its spread Ar has qualitatively 
similar consequences (see Fig. H), only that Rq at sufficiently large again reduces to the 
same small variance of the time distribution T(r). Note that, at small K^, Rq increases 
both with increasing At] and increasing Ar; this supports our claim that it is "sensitive" 
to the total longitudinal extension of the source Az ^ 2ro sinh At]. However, the relation is 
not linear (in particular in our numerical results a doubling of Ai] is seen to have less effect 



^Actually, the proper way of looking at this increase is by studying the variance of the function 
T(r): according to Eq. ( [All ) it increases from 0.89 fm/c for Ar = 1 fm/c to 1.49 fm/c for Ar = 2 



fm/c. 
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than a doubling of tq), making it hard to use in practice. It is obvious that the longitudinal 
velocity gradient (which decreases by a factor 2 when tq is doubled) has a stronger influence 
on Rq and R\\ than the geometrical width in t]. 



D. Transverse flow, M^-scaling, and kaon interferometry 

In this subsection we compare pion and kaon correlation functions. We discuss the M_l- 
scaling of the YKP radius parameters and its breaking by transverse collective flow. At 
the end of the subsection we formulate a program how to extract transverse flow from the 
M^-dependence of the YKP radius parameters. 

In Fig. ^we compare, for central rapidity pairs Y^^^^ = 0, the three YKP radius parameters 
for pion and kaon pairs, as functions of K_i_. The left column shows a source without 
transverse expansion, in the right column the transverse flow rapidity was set to rjf = 0.6. 
The onset of transverse flow has two qualitative effects: (i) the transverse radius acquires 



a i^_L-dependence 0, and (ii) Rq and y (t'^) begin to deviate from each other, as discussed 
in the previous subsection. (The equality i?y = {z^) remains exact because we are studying 
pion pairs at Y,-.^ =0.) The effects of flow on and (?) are seen to be weak, for both 
pions and kaons. 

Note also that at small K± the kaon radii are generically smaller than the pion radii, 



with or without transverse flow. This is also seen in experiment p6| , |38| . However, except for 
the change in the rest mass we have changed no parameters in the emission function, so the 
difference must be entirely kinematic. Indeed, it just reflects the fact that for thermalized 



sources like ( |3.1| ) the leading dependence on the particle rest mass is through the variable 



Mj_ = ym^ + K'j_. As discussed after Eq. ( |3.5| ) the source (|3.1|) depends only on M± if it 
does not expand transversally {i]t{r) = 0). In this case the correction terms in Eqs. ( p.l9| b,c) 
vanish exactly, and R\ = (y^), Rf = (5^), and -Rq = (t^) are also functions of M± only. This 



is shown in the left column of Fig. 10 where the YKP radii for a source without transverse 



flow are plotted as functions of M± and seen to exactly coincide for pions and kaons in the 
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common M^-range. Since in the absence of transverse flow the Boltzmann factor in ( p.l|) 
has no transverse gradients (we have assumed a constant temperature T), the transverse 
radius R± is M± independent and equal to the transverse geometric (Gaussian) radius R. 
It was pointed out in Refs. that transverse temperature gradients can also cause an 

M^-dependence of the transverse radius R±; but since the source remains in this case a 
function of M± only, the M_L-scaling of the YKP radii persists; it can only be broken by 
transverse flow. 



The breaking of the M^-scaling by transverse flow is shown in the right column of Fig. |Ty, 
for rjf = 0.6. It has two origins: the emission function itself is no longer a function of M± 
only (see (|3.5|) ), and the now non- vanishing correction terms in ( p.l9| b,c) depend on f3 and 
thus on both M_|_ and the rest mass m. It is obvious that the scaling violations induced by the 
pion-kaon mass difference are weak and require very accurate measurements. Furthermore, 
one may be worried that resonance decay contributions to the correlation radii [jl4[ (which 
we haven't discussed here) lead also to a breaking of the Mj_ scaling, because they affect 
pions more than kaons, and this may make it difficult to isolate the transverse flow effects. 
We refer to the detailed discussion of resonance decays in the context of the model source 
(p.lD in Ref. [|T7|. That study shows, however, that their influence on the M_|_-dependence 



of the transverse radius parameter i?_L is weak ||16| , p!7| . Furthermore, resonances tend to 
increase all three HBT radii (in particular the effective lifetime Rq), while the M^-scaling 
violations from transverse flow have the opposite sign for R\\ and R±, Rq. 

Detailed dynamical studies of the freeze-out process have shown that the transverse 
gradients of the temperature across the freeze-out surface tend to be small [pOlEOi. So the 



experimentally observed M^-dependence of the transverse radius p5|-[27|l is presumably due 
to transverse flow |4T|]. It was shown in Ref. that the strength of the M_|_-dependence of 
increases monotonously with the strength rjf of the transverse expansion. Alber has 
suggested to quantify the strength of collective flow by fltting the HBT radii to a power law 
in Ml , 
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RAM^) oc M^"^ , Ri{M^) oc MJ"' , (4.2) 

and using the magnitude of the extracted (negative) power as a flow measure. He found 
ai ~ 0.5 for Ri and smaller values for a±, with a tendency to increase for larger collision 



systems [41]. He interpreted this as a signature for strong longitudinal and weaker transverse 



expansion, the latter becoming more important for larger systems. 



In Fig. |Tl| we study the possible conclusions from such an exercise when applied to the 
results from our model. The left column shows double logarithmic plots for R± and R\\ as 
functions of M^. Obviously the assumption of a power law dependence is well justified for 
R\\ but somewhat marginal for R±. Rq{M^) cannot be approximated by a power law at all. 
In the right column we show the extracted powers as a function of ?7/, the scale parameter 
for the transverse flow. Since R± is not well represented by a power law, the extracted slope 
depends somewhat on the fit region, as indicated for the two sets of curves in Fig. |lT]b. 
Altogether it is, however, clear that for pions the power a± increases approximately linearly 
with T]f and for kaons somewhat more strongly. But even for large transverse flow rapidities 
rjf ~ 0.5 the power remains below 0.2. In contrast, the corresponding power ay in a fit 
R\\{M±) oc M_|_°" is already 0.55 in the absence of transverse flow, reflecting the strong 
boost-invariant longitudinal expansion. (Note that the decrease of -Ry with increasing M± 



is faster than the yT/Mj_-law suggested in Ref. - see also Ref. PJi^ for a discussion 
of this point.) As the transverse flow is switched on, ay changes much more weakly than 
a±, showing that R\\ is mostly sensitive to the longitudinal flow while i?_L is only affected by 
transverse expansion. Again, kaons are affected by the transverse flow more strongly than 
pions. 



V. CONCLUSIONS 

We have presented a numerical study of the Yano-Koonin-Podgoretskii fit parameters for 
the two-particle correlation function. Our starting point were the recently derived model- 
independent expressions for the YKP parameters in terms of second order space-time vari- 
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ances of the source emission function. These expressions allow for an easy evaluation of the 
YKP parameters as functions of the pair momentum K and for detailed parameter studies. 
We exploited them for a class of hydrodynamic models describing locally thermalized and 
collectively expanding sources, and we studied the dependence of the YKP parameters on 
the longitudinal and temporal extension of the source as well as on its longitudinal and 
transverse expansion velocity. 

In the context of such models it has been argued previously that in a certain approxi- 
mation (which becomes exact in the absence of transverse x-p-correlations of the source as, 
e.g., induced by transverse expansion flow) the YKP parametrisation achieves a perfect fac- 
torisation of the longitudinal and transverse spatial and the temporal extensions ("lengths 
of homogeneity") of the source, in the comoving frame of the emitting fluid element. The 
velocity of this emitting fluid element is then given by the fourth YKP parameter (the YK 
velocity). Here we have shown numerically that, within these models, these features are 
preserved even in the presence of transverse flow. The transverse radius parameter Rj_(K) 
gives the effective transverse size of the source, the longitudinal radius parameter R\\ (K) its 
effective longitudinal size, both in the local rest frame of the emitter as seen by pairs with 
momentum K. Also, the YKP parameter Ro(K) provides a direct estimate of the effective 
emission duration for particles with momentum K; this estimate is quite accurate as long 
as the average transverse flow rapidity remains below 0.5 and the "lifetime parameter" At 
is not too small (Ar > 1 fm/c). 

We also showed analytically and numerically that the YK velocity obtained from the 
YKP fit is indeed approximately equal to the longitudinal velocity of the emitting fluid ele- 
ment (the LSPS velocity), and that the small differences between these two velocities can be 
understood quantitatively in terms of asymmetries of the source around the point of maxi- 
mum emissivity (its "saddle point"). This enables us to interpret the rise of the YK rapidity 
with the pair rapidity F as a direct consequence of the longitudinal expansion of our source, 
and the M^-dependence of the slope of the function Y^^ (Y) as a signature for the thermal 
smearing of the particle momenta in the fluid rest frame. In Ref. we further showed that 
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the slope of the function Y^j^(Y) is nearly independent of the transverse flow rapidity rjf. 
Thus the YK rapidity is manifestly dominated by the longitudinal expansion and hardly 
affected by the transverse expansion at all. The latter causes, on the other hand, an M±- 
dependence of the transverse radius parameter R±, which in turn is completely unaffected 
by the longitudinal expansion. In the present model study, the YKP parametrisation thus 
not only leads to a factorisation of the (transverse and longitudinal) spatial and temporal 
aspects of the source geometry, but it also cleanly separates its transverse and longitudinal 
dynamics. 

The sensitivity of the YKP radius parameters to the transverse expansion of the source 
was investigated quantitatively in Sec. |1V D| . While the longitudinal radius parameter R\\ is 



affected very little by the transverse flow (its strong M_|_-dependence arises from the strong 
longitudinal flow), the transverse radius shows a considerable dependence on rjf (but none 
to the longitudinal flow). Furthermore, transverse flow breaks the exact M_L-scaling of the 
YKP radius parameters which we showed to exist for rjf = (see also [^). As explained in 



Sec. |1V L)| , both effects can be combined for a quantitative extraction of the mean transverse 



expansion velocity from the YKP radius parameters. The results from a comprehensive 
analysis of resonance decay contributions to the correlation function indicate that they 



don't jeopardise such a program. 

We would not like to close without remarking that the class of source models studied here 
is restricted in one crucial aspect: if the source is "opaque" |jl3[, (i.e. the particle emission 
is strongly surface dominated rather than being distributed over the whole emission region 
with a Gaussian geometric weight as assumed here), the resulting differences between the 
variances (x^) and (y^) lead to much larger contributions in Eqs. ( p.lSj) , ( p.l9|) and ( |4.1j ) than 
found in the present study. In Ref. |T5| it is shown that this affects strongly the interpretation 



of measured YKP parameters, in particular of the "temporal" parameter R^, which typically 
becomes strongly negative for small K± due to the negative contribution from (5:^ — y'^) in 
this case. However, that study also shows that presently available data on YKP radii from 



Pb+Pb collisions do not show any such indications for opaqueness of the source and 
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favor models with volume dominated emission. The NA44 data for HBT radii from 
heavy ion collisions with Pb targets at the CERN SPS, which within the standard Cartesian 
parametrization seem to be consistent with Ro = Rs, are more difficult to interpret because 
of the non-trivial shape of the acceptance window of this experiment in the Y — K± plane 
and the lack of information on the i^-dependence of these parameters. Also, the important 
consistency relations (|2.16|) so far have only been checked by the NA49 collaboration ||27 



whose data are in qualitative agreement with the numerical results presented here [^]. A 
quantitative comparison with the experiments will be presented as soon as finalized data 
become available. 
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APPENDIX A: SPACE-TIME MOMENTS O THE EMISSION FUNCTION 

Using cylindrical coordinates r, 77, r, with d^x = t dr drjr dr d(f), we can write the emis- 
sion function (B.ll) as 



S{x,K) d^x = T{t) P{r, 0) H{r, r]) dr dr] dr dcj) , (Al) 

with 



T{t) = , exp — ^ , (A2) 



P(r,0) = — e^('^)™^'^, (A3) 
27r 

H{r^ v) = j£y, exp (-^ - ^|^) cosh(r, - Y) e-«™^'^(-^) , (A4) 
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where we defined 



a{r) = ^cosh77t(r) , 
b{r) = ^smhrit{r) . 



The and r integrations can be done analyticaUy. We use 

d4) 



and define 



Jo ZTT 

/oo 
T(t) dT = To, 
-oo 

/oo 
rT(r)dr = ro^ + (Arf , 
-oo 

/oo 
T2r(T)c^T = To^ + 3To(AT)^ 
-oo 



The variance of the r-distribution T(t) is 

V To / 



Defining further 

..... _ r dr d77 g(r, 77) /o(&(r)) /(r, 77) 
\/ ir, v))^ ^^^^ ^^^^^^^^ 

we find for the non- vanishing moments up to second order 



2\ V mr)), 



-to 
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T 



{xt) 
(zt) 
(xz) 



Til , h{h{r)y 
' r cosh 77 



(smhrjcoshr))^ , 
h{b{r)y 



To 
T2 
To 

— { r sinh ri 
To\ 



(A19) 
(A20) 
(A21) 
(A22) 



For r]f ^ the (. . .)*-averages have to be done numerically. For convenience we also give 



{xi) 



T\ 
To 



r cosh r] 



h{h{r)r 



mr))L \ 
2 hm) \ i jMr)y 
mr))l. 



mr)). 



(cosh 7;)^ 



(A23) 
(A24) 



\ mr))/. 

Please note that for small arguments In{b) ~ 6". Thus for small r]f and/or K±, (xi) and 
(x^ — y^) vanish linearly and quadratically, respectively. 
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FIGURES 

FIG. 1. (a) The rapidity of the YK frame as a function of the pair rapidity Y (both measured 
in the source CMS), for pions (sohd) and kaons (dashed) and for transverse momenta K± = 1 MeV 
and K± = 1000 MeV. The transverse flow was set to 7]f = 0.3. (b) Same as (a), but shown as a 
function of M± for different values of Y. (c) The difference Yyk ~ ^lsps (^^^ text), plotted in the 
same way as (a), (d) Same as (c), but shown as a function of M_|_ for different values of Y. 

FIG. 2. The correction terms (x^ — y^) (a) and (xi) (c), for pairs with y = in the CMS, as 
functions of K± for different values of r]f. The third correction term (xz) vanishes at y = 0. Solid 
(dashed) lines refer to pion (kaon) pairs. Figures (b) and (d) show the same quantities, but scaled 
by the appropriate inverse powers of f3± (see text). 

FIG. 3. Same as Fig. but now for a fixed transverse flow r]f = 0.3 and pions only, but for 
two different pair rapidities, Y,^^ = (solid) and Y^^^ = 3 (dashed). The curves were evaluated in 
the YK frame. 

FIG. 4. Dependence of the correction term {xi) {K±) on the longitudinal width parameter Arj, 
for pion pairs with Y,^^ = 0. 

FIG. 5. (a) Rq and ^ Ij?) as functions of K± for rjf = and rjf = 0.6, for pion pairs with 
rapidity Y,^^ = and Y^^ = 3. The lifetime ^ (P) is evaluated in the YK rest frame, (b) Same 
as (a), but for i^ii and the longitudinal length of homogeneity y/ {z^) in the YK rest frame. For 



■^CM — ^11 ^'^'^ V {z^) agree exactly because /?/ = in the YK frame. 

FIG. 6. (a) Rq and \J (t^) as functions of Kj_, for different longitudinal gaussian widths Ajy. 
The diagrams are for pion pairs with Y^^^ = and for transverse flow rjf = 0.3. The lifetime (?) 
is evaluated in the YK rest frame which coincides here with the CMS and the LCMS. (b) Same as 
(a), but for Since Y^^ = 0, and the longitudinal length of homogeneity in the YK frame 
y/ {z^) agree exactly. 
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FIG. 7. Same as Fig. ^, but for different widths Ar of the proper time distribution in the 
source. 



FIG. 8. Same as Fig. ^, but for different average freeze-out times tq. 

FIG. 9. The YKP radii R± (top row), iiy (middle row), and Rq (bottom row), for Y^^ = 
pion (sohd) and kaon (dashed) pairs, as functions of K±. Left column: no transverse flow. Right 
column: transverse flow ijj = 0.6. In the lower right panel we also show the effective lifetime in 
the YK frame (P) for comparison. For more discussion see text. 

FIG. 10. Same as Fig. ^ but plotted as functions of M±. 

FIG. 11. (a) R± as a function of M± at Y^^ = 0, for pions (solid) and kaons (dashed) and 
different transverse flow rapidities rjf. (b) The scaling coefficient a± defined by R± « MJ"^ for 
pions (solid) and kaons (dashed) as a function of the transverse fiow rapidity r]f. The different 
results obtained by fitting in the regions M_i_ — m < 500 MeV/c^ and M± — m < 1000 MeV/c^ are 
shown separately, (c) Same as (a), but for (d) Same as (b), but for ay. 
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